Bulk Topological Invariants in Noninteracting Point Group Symmetric Insulators 
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We survey various quantized bulk physical observables in two- and three-dimensional topological 
band insulators invariant under translational symmetry and crystallographic point group symme- 
tries (PGS). In two-dimensional insulators, we show that: (i) the Chern number of a C-invariant 
insulator can be determined, up to a multiple of n, by evaluating the eigenvalues of symmetry 
operators at high-symmetry points in the Brillouin zone; (ii) the Chern number of a C„-invariant 
insulator is also determined, up to a multiple of n, by the C„ eigenvalue of the Slater determinant of 
a noninteracting many-body system and (iii) the Chern number vanishes in insulators with dihedral 
point groups D„, and the quantized electric polarization is a topological invariant for these insula- 
tors. In three-dimensional insulators, we show that: (i) only insulators with point groups Crt , C n h 
and S n PGS can have nonzero 3D quantum Hall coefficient and (ii) only insulators with improper 
rotation symmetries can have quantized magnetoelectric polarization P3 in the term P3E • B, the 
axion term in the electrodynamics of the insulator (medium). 



The study of novel topological phases of matter has 
become one of the most active fields in condensed mat- 
ter physics 1-21 . These phases are interesting because 
while deviating qualitatively from the conventional insu- 
lating phase, they cannot be described by any local order 
parameter in the Ginzbcrg-Landau- Wilson spontaneous 
symmetry breaking paradigm. Hcuristically, the word 
'topological' implies the presence of some global prop- 
erty, i.e., contributed by all electrons in the system, that 
distinguishes this special phase. Such a property is usu- 
ally marked by a global observable that takes different 
values in a topological phase and in the conventional in- 
sulating phase (normal phase) adiabatically continuablc 
to the atomic limit. In addition, 'topological' also implies 
that this global observable is quantized, or discretized, 
so that a topological phase cannot be adiabatically con- 
nected to the normal phase. Any quantized global ob- 
servable characteristic to a topological phase is called a 
bulk topological invariant. 

To date, topological phases have been identified in 
various systems, including those with intrinsic topologi- 
cal order (fractional quantum Hall states) 22-24 , topologi- 
cal band insulators 1-3,7 , topological superconductors 25-27 
and topological semimetals 18-20 . The focus of this paper 
is on topological band insulators, or simply topological 
insulators (TI). A typical, and first, example of TI is the 
integer quantum Hall state 28 . The bulk topological in- 
variant of this state is the quantized Hall conductance, 
o~xy = ne /h, where n £ Z, and its value remains fixed 
when magnetic field and gate voltage change within a 
certain range, resulting in a series of plateaus. Thou- 
less et al. 29 showed that the quantized Hall conductance 
is proportional to the Chern number (TKNN number) 
of a U(l) bundle over the 2D magnetic Brillouin zone 
(BZ), first time linking a quantized physical quantity to 
a topological number previously studied in the context of 
algebraic topology and differential geometry. It was then 



realized that since the Chern number is well-defined in 
any translationally invariant insulator regardless of the 
magnetic field, nonzero quantized Hall conductance can 
also appear without external magnetic field (or net mag- 
netic flux). Insulators with nonzero Chern numbers are 
later dubbed Chern insulators (or quantum anomalous 
Hall insulators) 30 . In general, the 2D Chern number of 
a translationally invariant insulator is given by 29 
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where |u n (k)) is the Bloch wavefunction of the nth band 
at k, and integral over BZ should be replaced with the 
one over magnetic BZ if there is an external magnetic 
field. From Eq.(l), it is clear that the Chern number is a 
global quantity contributed by all occupied electrons in 
the system. 

If on top of translational symmetry, there is time- 
reversal symmetry (TRS) in the system, then one can 
prove that the Chern number defined in Eq.(l) always 
vanishes, a case often referred to as trivial. In the absence 
of nonzero Chern number, the natural question then is: 
what is the bulk topological invariant for this system, if 
any? The answer is obvious when the z-componcnt of the 
total spin, S^ ot is conserved (commutes with the Hamil- 
tonian), because in this case the two spin components 
are decoupled so that we can define the Chern numbers 
for spin up (C-f) and spin down (Cj.) separately. Due to 
TRS, = — C±, hence the quantity 
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is a nontrivial bulk topological invariant, also known as 
the spin Hall conductance 2 ' 3 (when multiplied by K). 
However, when Sf ot is not conserved, the definition of the 
bulk topological invariant is not a simple extension of the 
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Chcrn number. According to Fu and Kane 7 , the nontriv- 
ial topology in a 2D insulator with translational symme- 
try and TRS can be characterized by a ^-number, i.e., a 
quantity that takes only two values. The explicit expres- 
sion of this i?2-number is not given here but we point out 
that it is linked to a quantized physical quantity, namely, 
the change in the 'time-reversal polarization' over half of 
the BZ. The -Z2 bulk topological invariant thus defined 
can be generalized to 3D insulators with TRS 7,16 ' 31,32 , 
where it is given by 

P 3 = / cPkeijkTriiFij - ^iA l A J )A k ] mod 1.(3) 

In Eq.(3), we have defined the non-Abelian Berry con- 
nection 

A? n (k) = i(u m (k)\di\u n (k)), (4) 

where |w„(k)) is the Bloch state and the non-Abelian 
Berry curvature or field strength 

J^.(k) = diAj - djAi + i[Ai, A,}. (5) 

From Eq.(3), we can see that P 3 is a global quantity 
depending on the band structure of all occupied bands in 
the whole BZ. In 3D, P3 has another physical meaning: it 
is the coefficient of the axion term of the electrodynamics 
of the insulator in the compact space (having periodic 
boundary conditions in all directions), 

H axion = P3E • B, (6) 

and is in principle measurable with electromagnetic in- 
duction experiments. In order for P3 to be a ^2-indcx 
(taking only two possible values), it must be quantized 
by some symmetry. In an insulator with no symmetry, 
P3 defined in Eq.(3) is an arbitrary number between zero 
and unity; but TRS quantizes it to either zero or one-half. 

From these previous studies, we can see that symme- 
tries of a system play a key role in the definition of its 
bulk topological invariants. TRS has a two-fold effect: it 
on one hand constrains the Chern number to be zero, and 
on the other hand makes possible the definition of a new 
Z2 number. An interesting question then is: if we con- 
sider adding other symmetries to an insulator, will there 
be new bulk topological invariant (s) characteristic to the 
added symmetries? Answering this question on a general 
ground is not easy, because it takes a detailed study of 
the homotopy groups of continuous mappings from the 
d-dimcnsional BZ to a target space of projectors under 
the symmetry constraints. Along this line, Schnyder et 
al. gave a list of all bulk topological invariants for each of 
the ten classes of Hamiltonians in the presence / absence 
of TRS, particle- hole symmetry and chiral symmetry for 
d= 1,2,3 27 . 

All the three symmetries considered in Ref.27 are sym- 
metries on internal degrees of freedom, which do not 
change the position or orbital character of an electron. 
Topological properties protected by these internal sym- 
metries are in general robust against random factors such 



as disorder. Nevertheless, the requirement of TRS ex- 
cludes a spectrum of interesting materials that have mag- 
netism while particle-hole symmetry and chiral symme- 
try are not to be found in realistic normal insulators but 
only in superconductors or maybe optical lattices. There- 
fore, it is desired to extend the study of topological in- 
variants to systems with other symmetries. In condensed 
matter electronic systems, it is natural to consider the 
lattice symmetries in the classification of insulators, be- 
cause of their universal existence in solids, nonmagnetic 
and magnetic. 

In 3D there are 230 types of lattice symmetries, in 
which the simplest, yet nontrivial, one is the inversion 
symmetry, which sends an electron at r to — r leaving 
its spin invariant. Insulators with inversion symmetry 
have been previously studied 33-35 and it is shown that 
this symmetry affects the topological properties in two 
ways. First, its existence places constraints on the values 
of bulk topological invariants such as the Chern number 
and the Z2 number: in Ref.34 it is shown that the par- 
ity (even or odd) of the Chern number is determined by 
the inversion eigenvalues of occupied bands at four inver- 
sion invariant k-points in the BZ, regardless of the Bloch 
wavefunction at all other k-points; the Zi number of a 
TI with TRS can also be determined by the inversion 
eigenvalues at all inversion invariant k-points 33 . Second, 
inversion symmetry alone brings about new topological 
invariants without TRS. We have mentioned that P3 is 
quantized to zero and one-half by TRS, and Ref.34 and 
35 shows that in the absence of TRS, P3 is still quantized 
by inversion symmetry and is equal to, up to an integer, 
half of the winding number of the inversion sewing matrix 
(to be defined in Sec.I(B)). 

In this paper, we study topological insulators subject 
to a more general class of lattice symmetries: crystallo- 
graphic point group symmetries, or simply point group 
symmetries (PCS) 36 ' 37 . Insulators with PCS are in- 
variant under a certain set of rotations and reflections 
that leave at least one point fixed in space. In the lan- 
guage of crystallography, by only considering PCS and 
lattice translational symmetry, we restrict the discus- 
sion to insulators with symmorphic space groups. Sim- 
ply speaking, in a symorphic space group, there are only 
point group operations plus translations by lattice vec- 
tors, while the a general space group can contain com- 
bined operations of a point group operation followed by 
a translation of fractions of lattice vectors. We ask two 
questions regarding the relation between PGS and bulk 
topological invariants: (i) how does a given PGS con- 
strain the value of a topological invariant such as the 
Chcrn number? (ii) can a given PGS give us new bulk 
topological invariants, or quantized global quantities? 

For 2D insulators, we answer these questions by an 
exhaustive discussion of all nine nontrivial PGS invari- 
ant insulators, divided into two types: those invariant 
under cyclic PGS C„ = 2.3,4.6 and those under dihedral 
PGS P^ n= i.2,3.4,6- For C„-invariant insulators, wc first 
show that the Chern number is determined by, or con- 
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strained by, the C m eigenvalues of each occupied single 
particle band at discrete high-symmetry points in BZ, 
where m evenly divides n. From there we further prove 
that the Chern number modulo n is determined by the 
C n eigenvalue of the many-body ground state, or Slater 
determinant in the noninteracting case. For D„-invariant 
insulators, we prove that while the Chern number is 
constrained to zero, the electric polarization, or more 
strictly, the dislocation of the electronic charge center 
from the nearest lattice point, becomes a bulk topological 
invariant, the value of which is determined by eigenvalues 
of C m . 

In 3D, there are in total 31 nontrivial PGS, and in- 
stead of exercising an exhaustive study, we focus on es- 
tablishing a link between the given PGS and the pres- 
ence/absence of two experimentally interesting topolog- 
ical invariants: nonzero quantized 3D Hall conductance 
and quantized magnctoelectric coefficient P3. We find 
that only C„, C n %, S n PGS are compatible with nonzero 
Hall conductance while all other PGS lead to zero Hall 
conductance in all three components. Moreover, we prove 
that the sufficient and necessary condition for the quan- 
tization of P3 (to zero and one-half) is that the PGS 
contains at least one improper rotation, and 2P3 can be 
expressed as the winding number of the sewing matrix 
associated with that improper rotation symmetry. This 
excludes the following point groups from having quan- 
tized P3: C„, D n , O, T, 11 distinct point groups in total. 

The paper is organized as follows. In Sec. I, we in- 
troduce concepts and basic formulae constantly used 
throughout the paper, including a brief introduction to 
Hilbert space representations of crystallographic point 
groups in Sec. 1(A), definition and basic properties of the 
sewing matrix associated with a PGS operator in Sec. 1(B) 
and the definition of basic properties of path-ordered loop 
integrals in BZ in Sec. 1(C). In Sec. II, we show how the 
Chern numbers in C n -invariant 2D insulators can be de- 
termined by eigenvalues of symmetry operators at high- 
symmetry k-points in Sec. 11(A), and how the Chern num- 
ber in C„-invariant 2D insulators can be determined by 
the C n eigenvalue of the Slater determinant in Sec. 11(B). 
In Sec. Ill, we show that D„-invariant insulators have zero 
Chern number in Sec.III(A) and that the electric polar- 
izations are quantized thus being topological invariants 
and how to calculate them in Sec.III(B). We study 3D 
insulators with PGS in Sec. IV, including the 3D quan- 
tum Hall effects in Sec.IV(A) and quantization of 3D 
magnctoelectric coefficient P3 in Sec.IV(B). In Sec.V we 
conclude the work. 



I. PRELIMINARIES 

Before proceeding with the relevant background mate- 
rial, we briefly clarify the notations we use throughout 
the paper. When used as indices, the Greek letters will 
always span the orbitals, while Roman letters denote ei- 
ther the bands or spatial directions (x, y, z or k x , k y , k z ). 



Throughout the text there are three different linear 
spaces in which operators are considered. First is the 
single-particle Hilbert space. The dimension of the single 
particle Hilbert space is N x N or b, where TV is the total 
number of sites and N or i, total number of atomic orbitals, 
including spin, per site. The second linear space is the 
orbital space. The dimension of the orbital space is N or i,. 
The third linear space is the occupied-band space. The 
dimension of this subspace is iV occ , where N occ < N or t, is 
the number of occupied bands. 

Summation of repeated indices is implicitly indicated 
when and only when the summation is over an index 
denoting the Cartesian coordinates: x, y, z and k x , k y , k z . 

A. Single-valued and double-valued representations 
of a double point group 

In physics, properties of Hilbert space representations, 
i.e., operators, of symmetry group are more relevant than 
those of the group itself. For a point group, before defin- 
ing the operators corresponding to each symmetry oper- 
ation, it is necessary to extend the point group, G, to 
the double-point group G D defined as G D = GL)(E*G), 
where E is a rotation of angle 2ir about any axis in space, 
which satisfies E 2 — E, and Eg = gE for g e G. In 
this definition, we notice that a rotation of angle 2n is 
distinct from the identity element in G D . All represen- 
tations of G D can be classified into single-valued repre- 
sentations and double- valued representations, in which 
D{E) = D(E) and D(E) = —D(E) respectively. In a 
physical system, the one-particle Hilbert space represen- 
tation of G D is a single-valued representation when the 
particle is spinlcss or of integer spin, and is a double- 
valued representation when the particle is of half odd 
integer spin. For example of the C n point group, if C n is 
the single particle operator for the n-fold rotation (and 
we will use a hat to denote operators in the single parti- 
cle Hilbert space) , then we have C™ = I if the particle is 
spinless or of integer spin, and C" = —I if the particle is 
of half odd integer spin, where I is the identity operator. 
To put the relations in a single equation: (7™ = (— \) F I 
where F is two times the total spin of the particle. 

B. Sewing matrix associated with a point group 
symmetry 

It is assumed that our system is defined on a lattice, 
has translational invariance and is noninteracting: 

H= £ h(k) (7) 

k£BZ 

= ^(k)4(k)c,(k), 

a,PMeBZ 

where a, f3 denote the orbitals within a unit cell. When 
wc say a Hamiltonian is invariant under point group G, 
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we mean given R £ G, there is 

Rh(k) = h{Rk)R, (8) 

where Rk is the transformed wavevector k under R. 
For example, if R is an n-fold rotation C n about z- 
axis, C n k = (k x cos(27r/n) — k y sin(27r/n), k x sm(2w/n) + 
k y cos(27r/n), k z ); if R is a mirror reflection about the 
xy-plane M xy , M xy k = (k x , k y , -k z ). 

At each k, the eigenstates of h(k) are called the Bloch 
states, whose annihilation operators 7„(k)'s satisfy 

[7„(k),ft(k)]=^(k)7 n (k) > (9) 

where n = l,...,N or f, is the band index. Then one can 
arrange j n in the ascending order of E n (k), such that 
for a fully gapped insulator, we have E n (k) < Ef for 
n < -ZV occ < N or b for all k 6 BZ. Combining Eq.(8) and 
Eq.(9), we have 

[R ln (k)R-\ h(Rk)} = E n (k). (10) 

This equation shows that i?7„(k)i? _1 is a Bloch state op- 
erator at Rk with the same energy. Therefore in general, 

R^k)^ 1 = J2 {B R (k)) mnlm {Rk), (11) 

mGocc 

where m £ occ. Considering degeneracies, (B R (k)) mn 
is not proportional to 5 mn . The N occ -by-N occ matrix 
B R (k) is called the sewing matrix associated with PGS 
operator R, and we will keep using curly letters to rep- 
resent matrices in the space of occupied bands. Defining 
\ip n (k)) = 7 ? l l (k)|0), where |0) is the vacuum state with 
no fermion, we can also put the sewing matrix definition 
in the form 

(B R (k)) mn = (^ TO (iZk)|% n (k)), (12) 

for m, n £ occ. (Throughout the paper, when it is 
a Greek letter inside |)/(|, the symbol denotes a col- 
umn/row vector in the single particle Hailbert space.) 

An alternative and useful expression of the sewing 
matrix can be found by noticing that each Bloch state 
\ip n (k)) is the direct product of a plane wave e lk r and 
an orbital part \u n (k)), and noticing that the rotation 
operator R acts on the two parts separately sending e lk r 
to e iBkr and |u„(k)) to R\u n ): 

(B R (k)) mn = (u m (Rk)\R\u n (k)}, (13) 

, where R and all capital letters with a tilde are matrices 
in the orbital space and when a Roman letter appears 
in |)/ (|, the symbol denotes a column/row vector in the 
orbital space. 

For a given PGS operator R, there are points in BZ 
at which ki nv = Rki nv . These points are called high- 
symmetry points, or points of symmetry. At ki nv , from 
Eq.(8), we have 

[iU(k im ,)] = 0. (14) 



Therefore, we can find a common set of eigenstates of 
both R and h(ki nv ). In that basis, the sewing matrix is 
diagonal 

(BR(k inv )) mn = R m (k)S mn , (15) 

where R m is the eigenvalue of R on the mth band. The 
determinant of the sewing matrix at ki nv can then be 
calculated as, using Eq.(15), 

det[B R (k inv )} = Yl R n (k inv ). (16) 

nGocc 

Since the determinant is independent of the choice of 
basis, dct[B R (ki nv )] is a gauge invariant quantity. 

C. Monodromy 

Path-ordered loop integrals and line integrals of ex- 
ponentiated Berry connection are extensively used in 
Sec. 11(A). Here we briefly develop basic properties of 
these integrals. Alternative formulas can be found in 
Ref.[40]. 

A path-ordered loop integral of Berry connection de- 
scribes a unitary evolution of the linear subspace spanned 
by the N occ eigenstates of H (k) around a closed circuit 
in the BZ. Mathematically, it is expressed in terms of the 
non-Abclian Berry connection 

W L = Pexp(i j A{k) ■ dk), (17) 

where L means a closed loop and P means 'path ordered'. 
Here we point out that in 2D, the determinant of this 
integral 

det(Wi) = dct(Pcxp(i <j> A{k) ■ dk)) (18) 

= det(cxp(i <j) A{k) ■ dk)) 

= cxp(iTr{<j, A(k) ■ dk)) 
= cxp(i0 B ), 

where <pB is the Berry phase associated with the loop. 

In Section 11(A), we will calculate the determinants of 
loop integrals that encircle a portion of BZ. To calculate 
them, it is convenient to define the Wilson line between 
the occupicd-subspaces at ki and k2, obtained by expo- 
nentiating the Berry connection: 

W kik2 = Pexp{i / A(k)dk). (19) 

In a finite system with discrete k, we have, equivalently, 

(^kik 2 )mri = (20) 

^( Um (k 1 )| Uo (k / 1 ))K(ki)|u 6 (k^))( U6 (k^)|..>„(k 2 )), 

a, 6. . . 
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where k'^kj... form a path connecting ki and k2. 

This connection is by definition an N occ x N occ ma- 
trix. We can also define an orbital space operator t/kik 2 
associated with this matrix: 

u klk2 = (^ik»)yM k i)>< u j( k 2)l (2i) 

= Pcxp(i / P k dP k • dk), 

where P k = Sigocc l u «( k ))( u '( k )l i s tne projector onto 
the occupied subspace at k. In a finite size system, al- 
ternatively, one has 

^ kl k 2 =p kl ( n ^.)i 3 k 2 . (22) 

3=1,2,- 

The advantage of using J7 kl k 2 is its gauge invariance as 
Pk only depends on the Hamiltonian. If the Hamiltonian 
has a symmetry R, it can be proven that 

RUk^R 1 = URk 1 Bk 2 > (23) 

although no simple relation exists for their corresponding 
matrices in the occupied band space W kl k 2 and ^Rki_Rk 2 - 
This property will be extensively used in the monodromy 
proof provided in Section 11(A). 
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FIG. 1. Schematic of the loop integrals used in the mon- 
odromy proof of the relationship between the Chern num- 
ber and the eigenvalues at high symmetry points within the 
Brillouin zone in C n invariant topological insulators with 
n = 2,4, 3, 6 in (a), (b),(c), and (d) respectively. 



II. TWO-DIMENSIONAL INSULATORS WITH 
CYCLIC POINT GROUP SYMMETRIES 

There are in total 10 point groups in 2D. One is the 
trivial group containing only the identity, and among the 
nontrivial 9 PGS, there are four cyclic PGS denoted by 
Cn=2.3,4.6- The symmetry group C n is generated by an 
n-fold rotation about the z-axis that is out-of-plane. In 
this section we will study C n -invariant insulators in two 
dimensions. 



A. Chern number in terms of rotation eigenvalues 
at high-symmetry points 

In the following, we endeavor to prove in a C n -invariant 
insulator, the Chern number modulo n is given by C m 
eigenvalues of all occupied bands at each C m -invariant 
k-point for each m dividing n. The case of n = 2 has 
been studied in Ref.34 and the result is 

(-1) C = II UT)UX)UY)UM), (24) 

where Q(k = T,X,Y,M) is the eigenvalue of the oper- 
ator C2 at ^-invariant k-point on the ith band. This 
equation shows that the parity (even or odd) of the 
Chern number can be determined by knowing the inver- 
sion eigenvalues at four high-symmetry points. 

Now we move beyond the simple and already discussed 
case of C2 symmetric systems to calculate the Chern 



number in C n -invariant insulators for n = 3, 4, 6. We will 
calculate the Chern number by relating it to the deter- 
minant of a special closed loop integral of exponentiated 
Berry connection defined in BZ. In the main text the 
detailed calculation is only shown for ^-invariant insu- 
lators; for Ca^-invariant insulators, the main conclusions 
are stated here while details are provided in Appendix B. 

For n = 4, the loop is chosen as A = TXMYT shown in 
Fig. 1(b), and the determinant of the path-ordered loop 
integral becomes 

dct[Pcxp( j> A(k) ■ dk)} = e i7rC/2 . (25) 

To prove this we note that due to C4, we have the fol- 
lowing symmetry in the field strength 

Tr[F xy (k x , k y )\ = Tr[F xy {-k y , k x )\ (26) 

(for a proof see Appendix E) and hence the integral of 
Berry phase inside A is one fourth of the Berry phase 
in the whole BZ (2wC). On the other hand, the loop 
integral is given by the connection matrices defined in 
Sec.I(C) 

Wx=U V xU xm UmyUyt- (27) 

With C4 symmetry at hand and utilizing the periodic- 
ity of the BZ, we notice that CJJtyC^ 1 = Urx and 
C± x UymC± = Umx- We may further simplify the deter- 
minant as (see Appendix A for more additional details) 
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det(W A ) = dct(B Ci {T)U V yB C2 {Y^UymBc, (M)U M yU yt ) (28) 
= dct(B Ci (T)B c 1 2 (Y)BcM)) AzWtyUymUmyUyt) 
= det(B Ci (T)Bc 1 2 (Y)B Ci (M)), 



where we have used Uyu = Umy- From Eq.(25,29) and 
Eq.(16), we obtain 

* c II Ur)UM)cHY), (29) 

nGocc. 

where £i(k) is the eigenvalue of C4 at the ^-invariant 
k-point on the ith band. Since C\ = (— 1)^, we have 
( 2 (k = X,Y) = (-1) F or 

C 1 (Y) = (-1) F £(Y). (30) 

Substituting Eq.(30) into Eq.(29), we obtain 

* C = II (-l) F Zi(m(MMY). (31) 

Eq.(31) tells us that given the eigenvalues of C4 at T, M 
and those of C 2 at Y, one can determine the Chern num- 
ber up to a multiple of four. It may be disconcerting 
that in Eq.(31) we only include the C 2 eigenvalues at Y 
but not at X, but in fact they are identical: in systems 
with C4 symmetry, £j(0,7r) = C^ 71 ", 0) as the two points 
are related by C4 rotation. 

For n = 3 and n = 6, the proof takes very similar steps 
which we leave to Appendix B, and only differs in that 
one chooses the loops as shown in Fig. 1(c) and Fig. 1(d), 
respectively. We simply quote the salient results here 

e ^c/3 = -Q (-ifOiWOiWOiiK'), (32) 

i£occ. 

e *rC/3 = -Q (- 1 )F rh ( T ) .(K) Ci ( M ), (33) 

iGocc. 

where 9i (k) and rji (k) represent the eigenvalues of C3 and 
C*6 at C3- and Cg-invariant k-points, respectively. 
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B. Chern number in terms of the C n eigenvalue of 
the Slater determinant 

We have shown how the Chern number is related to 
symmetry eigenvalues of the single particle states at high- 
symmetry points. Now we will see that another con- 
nection between the Chern number and C„-invariance 
presents itself as one considers the C n eigenvalue of the 
many-body ground state (under the assumption that this 
is non-degenerate), or the Slater determinant state of a 
C ra -invariant insulator: 

cxp{iC— ) = f n (N)p n , (34) 

n 

where / n (iV) is either +1 or —1 depending on the to- 
tal number of sites and p n is the eigenvalue of C n of 
the Slater determinant state |$o), i- e -> C n \$o) = p n \$o}- 
Eq.(34) can be proved by expressing p n in terms of the 
eigenvalues of the single particle states at high-symmetry 
points and then using Eq. (24,32,31,33) to relate p n to the 
Chern number for n = 2, 3, 4, 6. 

We begin by looking at the case of ^-invariant insula- 
tors. The Slater determinant of the ground state is given 
by 

i$o> = n n ^( k )i°>- ( 35 ) 

keBZ iSocc. 

We can pair k and — k and this results in only four k's 
that cannot be paired as they are inversion invariant. So 
the Slater determinant can be written as 

l^o) = II ll(rH(Yhj(X)^(M) (36) 

i£Locc. 

x II (7l(k) 7 l(-k))|0). 

keBZ/2 

Note that in writing down the above equation, and 
throughout the rest of the Section, we have implicitly 
assumed that all high-symmetry points exist in the fi- 
nite system, and for a discussion of all other cases, see 
Appendix C. Now consider the action of (7 2 on this state: 



I 

C 2 |$o)= II (C 27 l(r)C 2 - 1 )((7 27 l(r)C 2 - 1 )((7 27 |(A)C 2 - 1 )(C 27 l(A/)(7 2 - 1 ) (37) 

i£occ. 

IJ (^(k^^C-k)^)!^) 

keBZ/2 

= (-l) 7V - (Ar - 4)/2 det(B C2 (r))dct(S C2 (A))det(B C2 (r))det(S C2 (Af)) [] dct(£ C2 (k)£ C2 (-k))|$ }. 

ke-BZ/2 
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Using the sewing matrix property to see that Bc 2 (k) = Combining Eq(38) and Eq(24), we have 
(-l) ir S^ 1 (-k), one obtains from Eq.(37) 

C2 (-if = (_l)(^-l)JWV/2 p2 (39) 

for ^-invariant insulators. 

The case with n = 4 can be similarly studied. Here 
four generic k's can be grouped with the exception of the 
(f-i)jv (jv-4) r an d M points which do not pair with any other k as 

C 2 |$ ) = (-1) ° r J! Ci(r)Ci(-X")Ci(^)Ci(Af)|$o)-they are invariant under C 4 . Additionally, X and Y, pair 

ieoec. between themselves to make a group of two. The Slater 

(38)dcterminant is then given by 



1*0)= n^^Wf^ 1 )^ 7 )) II (7l(^Kh}(-K,k x hH-k*,-kyh!(k y> -k x ))}\o) (40) 

ieocc. ke-B2/4,k^ki„„ 

= {-l)"™"'* det(B C4 (T)) det(B C4 (M )) dct(B C4 {X)B Ci (Y)) 

Bc 4 ( k x,k y ) B c 4 (~ k y^ k ^) B cA~ k x, -ky)B Ci (ky,-k x ))\^ )- 

keBZ/4,k^ki„„ 



Since C4 is a four- fold rotation, two consecutive CVs 
equal a two-fold rotation, and that four consecutive CVs 
equal a complete rotation. In terms of sewing matrices, 
these simple facts are represented by 

BcAX)BcAY)=Bc 2 (X)i (41) 
£(k)£(C 4 k)S(C 4 2 k)£(Cfk) = 

By substituting these relations into Eq.(40), we obtain 

cm = (-i) N - l{F - 1)N/i+F] n tiFMM)Ci(X)\* ). 

iGocc. 

(42) 

Combining Eq.(42) and Eq.(31), we have 

i c = (_l)^-i)^-«^p 4 . (43) 

Finally, we discuss systems possessing C3 symmetry 
and Cg in a straightforward fashion. Here, with the ex- 
ception of the three k's that are invariant under C3: T, 
K, and K', whose corresponding Bloch state operators 
7j(kj m) ) cannot be grouped in threes in the Slater deter- 
minant for each band, all other 7j(k)'s (i.e., at other k's) 
can be grouped in threes. Following the same procedure, 
one finds 

C 3 |$o) = 1] ( 44 ) 

iGocc. 

JJ det(B(k)B(C 3 k)B(C 3 2 k))|$ ). 

keBZ/3M^k inv 

Since det( J B(k)B(C 3 k)B(C|k)) = (-1) F because three 
Ca's equal a complete rotation, we have 

c 3 |$ ) = (-i)™-^- 3 )/ 3 J] 0i(r)0i(iOW)l*o>. 

(45) 



Therefore the Chern number can be expressed in terms 
of /9 3 as 

exp(iC^) = (-l)™^/ 3 p 3 . (46) 

For C'e invariant insulators, similar steps (not shown 
here) lead to 

exp(ic|) = (-1)™-"/%. (47) 

Up to this point we have proved Eq.(34) for all C n - 
invariant insulators. Physically, it shows that the Chern 
number modulo n in a C„ invariant insulator is ex- 
actly the total angular momentum of the ground state. 
Although the Chern number was originally defined in 
terms of the Berry connection of single particle states, 
one can generalize its definition to an interacting sys- 
tem using flux insertion, also known as the distorted 
periodic boundary conditions, if the ground state is 
non-degenerate. It is then straightforward to see that 
in a weakly interacting system, the relation between 
the Chern number and the total angular momentum of 
ground state still holds because both numbers are quan- 
tized and therefore cannot change infinitesimally as in- 
teraction is adiabatically turned on. This relation largely 
simplifies the calculation of Chern number in a weakly in- 
teracting system, because only ground state information 
with normal boundary condition is needed to obtain the 
total angular momentum, or the C„ eigenvalue. 

III. TWO-DIMENSIONAL INSULATORS WITH 
DIHEDRAL POINT GROUPS 

Out of 9 nontrivial PCS in 2D, five are called dihcral 
PGS, denoted by £^=1,2,3,4,6- The structure of the di- 
hcral group is simple: D n = C n U (C' 2 * C n ), where C 2 
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is a two-fold rotation about some in-plane axis (referred 
to as x-axis in later text). For later discussion, we define 
the basis vectors in the real space and reciprocal space in 
four out of five 2D Braivis lattices in Fig. 2: rectangular, 
centered rectangular, square and triangular lattices. The 
parallelogram lattice is excluded from the discussion as 
it cannot be D„-invariant. 

In 2D, _D„-group is identical to C„„-group for spinless 
particles, where C nv is generated by C n and a mirror 
plane of xz, M xz . This is because both C 2 and M xz will 
send (x,y) to (x, —y). For spin-1/2 particles, the C 2 acts 
as a z in the spin space, while M xz acts as a y in the spin 
space. Therefore, in principle a -symmetric 2D sys- 
tem is not a _D„-symmctric system. However, through- 
out the section below, they are completely equivalent, 
because we only need a symmetry that sends a state at 
(k x , k v ) to one at (k x , —k y ) of the same energy and both 
C 2 and M xz have this property. Therefore, all conclu- 
sions obtained here apply to C nv -symmetric insulators 
without any adaptation. 



A. Vanishing Chern number 

A dihedral PGS has a strong constraint on the Hall 
conductance of a 2D insulator: 




FIG. 2. Primitive cells of four 2D Braivis lattices: (a) rect- 
angle lattice, (b) centered-rectangle lattice, (c) square lattice 
and (d) triangle lattice. The basis vectors of the real lattice, 
ai : 2, are plotted as blue arrows and the basis vectors of the 
reciprocal lattice, bi^, are plotted as red arrows. All possible 
C'2 axes are plotted as dashed lines in (a-c). In (d), dashed 
lines and dotted lines represent two sets of C 2 axes. D3 PGS 
includes either one of the two sets and Dq PGS contains both 
sets. 



B. Quantization of the electric polarization 



a xy = 0. 



(48) 



Since in a nonintcracting insulator the Hall conductance 
is proportional to the Chern number, Eq.(48) implies that 
the Chern number of an insulator with any dihedral PGS 
must be trivial. A heuristic understanding of Eq.(48) is 
that since a xy changes sign under C 2 , if the system is C' 2 
invariant, a xy can only be zero. This simple argument 
can be put in a rigorous form presented in Appendix D. 
An alternative proof of the vanishing Chern number is by 
using the general transform property of held strength T xy 
Eq.(83) proved in Appendix E, and substituting R = C 2 
to obtain 

•F zy(k x , ~ky) — —Bc^(k x ,ky)J- z y (k x ,k y )Bg,(k x ,k y ). 

(49) 

Substituting Eq.(49) into Eq.(l), we hnd 



C 



- dk x [ Tr(F xy (k))dk y + / Tr(T xy (k)dk y ) 

,7r J— it J-n JO 

Tr{JF xy{k x , k y ) + J~ xy (k X7 kyjjdky 

(50) 



h 1 dki 

0. 



Although Eq.(48) and Eq.(50) are equivalent to each 
other in a noninteracting insulator, in an interacting in- 
sulator with dihcral PGS, Eq.(50) loses its meaning due 
to the absence of Bloch states, but Eq.(48) still holds. 



As the Chern number always vanishes in these insu- 
lators, we need some other bulk topological invariant 
to identify nontrivial _D„-invariant 2D insulators. When 
Chern number vanishes, one can always find a continuous 
gauge, or a choice of continuous Bloch wavefunctions on 
the whole BZ with periodic boundary condition. With 
this gauge choice it can be proved that the electric polar- 
ization, or the center of charge, is unambiguously defined 
within a unit cell, and its position vector is given by 



(51) 



where ai ; 2 are the unit cell basis vectors of the lattice 
and are in general within the range [0, 1)., given by 
the integral of the Berry connection 39 : 

p. = J_ / dkt [ dk 2 Tr[Ai(hbi + k 2 b 2 )}, (52) 
2lT Jo Jo 

where ki is the component of k in its linear expansion 
in reciprocal lattice vectors b lj2 . Here we briefly discuss 
why the vanishing of Chern number is necessary for the 
definition of polarization, taking i = 1 for example. First 
notice that J dkiTr[A\ (fei, k 2 )\ = P\(k 2 ) is the ID po- 
larization at fixed k 2 . If the Chern number is nonzero, 
this number changes by exactly C from k 2 = to k 2 = 1 . 
Therefore the integral in Eq.(52) depends on the integra- 
tion range of k 2 . For example, the integral takes different 
values for integrating in [5, 1 + S) and in [0, 1). Since no 
physical quantity should depend on the choice of BZ, the 
polarization in this case is meaningless. 
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TABLE I. Quantized values of electric polarization (pi,p2) in D„=2,3,4,6-fnvariant 2D insulators. From the first row: the value 
of n in D n , the transform matrix of n-fold rotation in the basis of bi, 2 , the unit cell in the Braivis lattice with all C 2 -axes, 
symmetry constraints on the values of pi j2 by C n , additional constraints by C 2 , all nontrivial (pi,p2) pairs, the visualization 
of the charge center (large dotted circle) in the lattice defined by solid dots. 



Before presenting the formal results of quantized elec- 
tric polarization in _D„-invariant insulators, it is helpful 
to establish an intuitive understanding. First we remem- 
ber that the electric polarization is quantized to and 
1/2 (in unit eh where L is the length of the system) in a 
ID insulator with inversion symmetry. This result can be 
easily understood in a heuristic way: suppose the system 
has polarization p, but due to inversion, the polariza- 
tion is also —p. Then since p is only well defined up to 
an integer, the only possible values of p are p = and 
p = 1/2. Now given a Z^-invariant system with polariza- 
tion (pi,P2), then from C4-invariance it must be equiv- 
alent to (— P2,Pi), and since both p\ and P2 are defined 
up to an integer, we have equations p± = — p 2 + integer 
and P2 = Pi + integer, the only solutions to which are 
Pi = p2 = and pi = P2 = 1/2. One is reminded, how- 
ever, that although we only used C4 to prove the quan- 
tization of pi,2j we implicitly assumed that the Chern 
number is zero, without which the polarization is unde- 
fined. (This important fact is obscured in this heuristic 
picture.) Following the same approach, one can prove 
that for Da^-invariant insulators, (pi,p2) is also quan- 
tized and the quantized values can be similarly written 
down. Despite being physically clear, this simple argu- 
ment does not serve as a proof of the our statement, and, 
further more, it does not give the actual value of (pi,P2) 
for a given band structure. In what follows, we prove the 
quantization of electric polarization in D„-invariant insu- 
lators analytically and express the quantized polarization 
in terms of C n - eigenvalues at points of symmetry. 



quantizes p\ j2 - Using sewing matrix Bc n , we have 
Pl = 4z I dMfoTr<u(k)|(%-^-C' n |u(k)) 



2tt 



(53) 



± J d* 1 dA a rr<«(C n k)|Bj, n (k)^ : Bcr B (k)|«(C r n k)) 
^- f dMfc27>MC„k)|^| u (C„k)> 



i / Tr[B Cn {^) ] 4~ B o n (k)]&i*2 
J dki 



Using the unitarity of Bc„ (k) , the integrand in second 
term of the last line can be rewritten as 

Tr[B Cn (k)t 4r B cM] = ^ln[dct(B c „(k))]. (54) 

Substituting Eq.(54) to the second integral of last line 
of Eq.(53), we find (taking i = 1 for example while the 
same proceeds for i = 2) 

; dk 1 dk2Tr[B Cn (ki,k2^-^B Cn (ki,k2)} 



i 

27 



dk2{ J dk\ 



dln[dct(g c „(fci,fc 2 ))] 
dk\ 



}■ (55) 



The integral over k\ will give 2qiiri where qi is the wind- 
ing number of dct[Bc n ], which is a pure phase, at a fixed 
k 2 . However, we have chosen a smooth gauge so qi(k 2 ) 
must be a smooth function of , which, in the case of a 
discrete function, is a constant. That is to say, the wind- 
ing at each k 2 must be the same. Therefore, we have 



/ dk 1 dk 2 Tr[Bc n (h,k 2 y^ r Bc n (k 1> k 2 )} 
2tt J dki 



Now we prove that C n PGS, a subgroup of D n PGS, 



dk 2 qi{k2) 
qi 6 Z. 



(56) 
(57) 
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Use Eq.(56) in the last line of Eq.(53), we have 

Pl = h [ k ^ T r(u(C n k)\^-\u(C n li)) + gi(C n )<.58) 

The first term of the right hand side of Eq.(58) may 
be further simplified using 



(u(C n k)\ — \u(C n k)} 

(u(C n k)\ U (C n k + dkiC n bi)) - 1 



(59) 



dki 



<" ) Mk , )l47l^(k'))|k'=c„k, 



<dk 



where RS n > is the 2-by-2 matrix describing the rotation 
of the basis vectors bi t 2 under C n , summarized in the 
second row of Table I. Inserting Eq.(59) into Eq.(58) 
and after some rearrangement, we obtain 



Pi 



(60) 



Eq.(60) is the general form of symmetry constraints by 
C n -invariance on the electric polarization. 

In addition to the constraints by C' n symmetry, C 2 
symmetry also puts constraints on (pi,p 2 ). These con- 
straints are derived in a way very similar to the way 
sketched in Eq. (53-58). The result is given in the same 
form of Eq.(60): 



(61) 



Where R ^ 2 > is the transformation matrix of the basis bi.2 
under the two- fold rotation C 2 (also listed in Table. I). 
One can see in many cases, this additional symmetry 
does not give any additional constraints/quantization of 
Pi,2, but there are a few important exceptions. For n = 1, 
i.e., when there is no rotation symmetry about z — axis, 
the C 2 symmetry gives quantization to the polarization 
perpendicular to the rotation axis (defined as cc-axis). 
This is because for this system, at each k x , H(k x ,k y ) is 
a ID inversion symmetric insulator whose polarization is 
quantized, and due to continuity of the gauge, the po- 
larization is the same for all k x . Therefore, the total 
polarization is quantized. For Z^-msulators, if the C 2 
axis coincides with the bisect of ai and a 2 , i.e., in the 
case of centered-rectangle lattice, the C 2 rotation sym- 
metry gives p\ — p 2 G Z as the additional constraint. 
Another important case in which the C 2 gives additional 
constraint is in I?3-insulators, when the C' 2 axis coincides 
with ai (or sl 2 ). There the constraint by C 2 symmetry 
forbids any nontrivial (non- integer) values of p\ and p 2 to 
exist. Row 4&5 of Table I list all constraints on (pi,P2) 
for insulators with Z?„ = 2 i 3 ; 4 i 6 _ i nvar iance on four of the 
five 2D Braivis lattices. Notice that for D3-invariant 
insulators, there are two possible alignment of the C 2 
axes in the triangular lattice: cither along ai (and its C3 
equivalent axes) or along the bisect of ai and & 2 (and its 



C3 equivalent axes). Applying the constraints, one can 
easily write down all nontrivial values of (^1,^2) 7^ 0, 
listed in the last row of Table I. Below each pair of 
(PiiP2)j the electronic charge center in the real space 
lattice is plotted. 

A few remarks are due regarding the results shown in 
Table I. First we emphasize that all results depend on 
the vanishing Chern number as a necessary condition, 
without which the 2D polarization is undefined. Second, 
although C 2 -invariance is a sufficient condition for the 
vanishing Chern number, it is not necessary for quanti- 
zation of (pi,p 2 ). In other words, given a 2D insulator 
with zero Chern number but without C 2 invariance, the 
constraints placed by C n shown in the fourth row of Ta- 
ble I arc still valid. Third, we notice that for Z^-invariant 
insulators with C 2 -axis aligned with ai.2 or Z^-invariant 
insulators, (pi,P2) = (0,0) by symmetry, i.e., there is no 
nontrivial polarization for these insulators. 

In Sec. II, we have shown the relation between the 
Chern number and symmetry eigenvalues at high- 
symmetry points in BZ. Can we also express the elec- 
tric polarization (pi,p2) in terms of the eigenvalues of 
symmetry operators in D n point group at high symme- 
try points? To have explicit values of pi >2 we need to 
evaluate the integer q on the right hand side of Eq.(60), 
defined through a sewing matrix in Eq.(56). From the 
definition, we can see that 2Trqi(C n ) is the phase differ- 
ence between the determinant of the sewing matrix at two 
k-points separated by b^. In principle, qi depends on i 
through b; and PGS. In practice, due to the presence of 
constraints, we will see that we only need to calculate 
qi(C 2 ) and (71(63) — q 2 (C3), while all other qr's can either 
be derived or are not needed. 

To calculate qi(C 2 ), first use its definition in Eq.(56) 
and the periodicity of the gauge to obtain 



Qi(C 2 ) 



2tt 



1/2 



dki 



1/2 



dln[det(B C2 (ki,0))] 
dki 



(62) 



In Eq.(62) we have taken k 2 = because the integral is k 2 
independent. Then we apply the sewing matrix property 
Eq.(F3) and have 



£ C2 (kH(-l) F 4 2 (-k), 



(63) 



which leads to 



ln[det(Bc 2 (fci,0))] = iN occ Fir - ]n.[det(B C2 (-h, 0))] 
dln[det(Bc 2 ( k u 0))] d\n[det(B C2 (~k 1: 0))] 



dki dki 
Substitute Eq.(64) into Eq.(62), and we have 

' 1/2 dln[det(S C2 (fci,0))] 



qi(C 2 ) 



dki ■ 



7T Jq dki 

-i , r dct(#c 2 pO) n 

T'"[d»IF |mod2 



(64) 



(65) 
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Eq.(65) can also be put in the form 

Wi(C 2 ) _ TT CnPO 



(-1)0 



n 



Cn(T) 



(66) 



The calculation of g2 follows exactly the same route only 
replacing X with Y. The result is 

Cn(Y) 



(_1)92(C 2 ) _ 



Cn(r) ' 



(67) 



We then use the constraints 2pi = qiiCz) and 2p 2 = 
1i (C2) for Da-invariant insulators to obtain the explicit 
values of P\$- 



(_!)*. = JJ 

(_i)^ = J| 



Cn(r) ' 
C™(r) ' 



(68) 



Mark that the inversion eigenvalues at M do not enter the 
formula due to the constraint of vanishing Chern number: 

n l60cc c.(r)opoc.onc*(A/) = i. 

9i(C2) can also be used in calculating pi in D 4 - 
invariant insulators, because D4 already implies D 2 . In a 
/^-invariant insulator, using the constraint 2p\ = 2pi = 
9i(C2) we hud 



( _ 1) 2 Pll2 = -Q 



n^occ. 



Cn(X) 

Cn(r) • 



(69) 



To calculate q\ (C3), we start from inspecting the phase 
of the determinant of the sewing matrix associated with 
C 3 : 



</>(k) = -ilndet(Bc 3 )(k) 



(70) 



in the BZ. Due to the continuous and periodic gauge 
choice, 0(k) is a continuous function satisfying 



(k + bO =0(k)+2 gi (C 3 )7r. 



(71) 



<t>{K) and 4>(K') can be determined by the C 3 eigenvalues 
at these two points 



<KK) = -iMU 0«(K)], 

n£occ 

4>(K') = -i\n{ l[ 6 n (K')}. 

n£occ 



(72) 



We can then express </>(k) at all six corners of the BZ in 
terms of <j)(K,K') and q\^{C 3 ) as shown in Fig. 3. Then 
we use the sewing matrix property 



£c 3 (k) = (-l) F 4 3 (C 3 k)^(C 3 - 1 k) 



(73) 



to have 



4>{K') + 2q 2 7t - 4>(K) = f Tr(B C3 (k)dBc 3 (k)) • dk 
(f dk+ f dk) ■ Tr(Bc 3 (k)dB C3 (k)) 

J A2 " A3 

2qnr - 2{<P{K') - <j>{K)), (74) 



where Ai 2,3 are integration paths marked in Fig. 3. From 
Eq.(74) we extract 

qi(C 3 ) - q2(C 3 ) = ^W<') - <KK))- (75) 

Eq.(75) and the constraints placed by C*3-invariance give 
the final expression of PiX- 



exp(z27rpi, 2 ) = 



n£occ 



e n {K) 

UK')' 



(76) 



4>k + 1(q 1 +q 2 )ii 




FIG. 3. The phase of det(Bc 3 ) at six corners of the BZ. Ai,2,3 
denote the paths of integration considered in the text. 

Up to this point, we have formulated the theory of 
quantized polarization in _D„-invariant insulators. In 
fact, we emphasize that the theory applies to a wider 
class of insulators: C„-invariant insulators with vanish- 
ing Chern number. The electric polarization serves as a 
complimentary quantum number when the Chern num- 
ber is zero (for example, in _D„-invariant insulators), and 
just like the Chern number, the electric polarization can 
be expressed in terms of the C„-eigenvalues of bands at 
points of symmetry. 



IV. POINT GROUP SYMMETRIC 
TOPOLOGICAL INSULATORS IN 3D 

In 3D, there exist 32 different point groups, each of 
which is generated by four types of operators: n-fold ro- 
tations C n , mirror plane reflections M and and rotation- 
reflection S n . S n is a combination of an n-fold rota- 
tion followed by a mirror reflection about the plane per- 
pendicular to the rotation axis, or S n = C n * M xy . In 
the following, we explore two fundamental responses of 
point group symmetric topological insulators: the 3D 
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anisotropic quantum Hall conductance and the magne- 
toelectric polarization. 



A. 3D Anisotropic Quantum Hall State 

The simplest 3D point groups arc the cyclic groups C n 
as they are comprised of only one rotation axis, which we 
assume to be the z-axis. The 3D single particle Hamil- 
tonian H(k x , k y ,k z ) can be constructed from a series of 
2D Hamiltonians Hk a ,k ;k. where k z is a momentum pa- 
rameter spanning — 7r to it. For every one of these in- 
dividual 2D systems, one can define a Chcrn number 
C{k z ). For an insulator, there is C(k z ) — constant, and 
if C{k z ) changes at some k z , there must be, generically, 
a 3D Weyl node in the bulk 18 " 20 - 34 . For n = 2,3,4,6, 
C(k z ) can be calculated using Eq.(24), Eq.(32), Eq.(31), 
and Eq.(33) respectively. Each 2D subsystem contributes 



a xy (k z ) = Ce ./h, and this allows us to easily calculate 
the total Hall conductance of the 3D system as 



T 3D 

'xy 



ch 



(77) 



where L z is the length of the sample in the out-of-plane 
direction and c is the lattice constant along the the same 
out-of-plane direction. On the other hand, the other 
two transverse conductances a xz and o~ yz must be zero 
as both a xz and a yz transform as vectors on sy-planc 
for a rotation about z-axis. Therefore they cannot take 
a non-zero value in a rotation invariant singlet ground 
state. Based on the preceding argument, as soon as there 
is more than one rotation axis, all three components of 
transverse conductance are quantized to zero. This di- 
rectly applies that D n , D nd , T hl T d , O and Oh results in 
Cxy = Cyz = o~ xz = 0, because they all have at least two 
rotation axes. 

Beyond rotational symmetries, we consider the effect 
of mirror plane symmetries. For C n h point group, where 
the mirror plane is perpendicular to the rotation axis (z- 
axis), we can still have a non-zero a xy given by Eq.(77). 
For C nv (n = 2,3,4,6, while C\ v = C\h) point group, 
since the mirror plane is parallel to z-axis, a xy = 0, be- 
cause o~ xy changes sign under the mirror symmetry while 
the ground state is a singlet eigenstate invariant under 
any operation in the point group. As a result, all com- 
ponents of Hall conductance vanish for C n „-invariant in- 
sulators. 

Finally we discuss the Spiegel group S n {n = 2,4,6, 
while S*3 equals C^h and has been discussed). 52 is the 
same as inversion, and we borrow the result from Ref. 34 
which states that all three transverse conductances may 
be quantized to nonzero values, depending on the Chern 
numbers on three planes defined by ki =x>y>z = respec- 
tively. S4 and Sq has C2 and C3 as subgroups, respec- 
tively, and hence have <j xz = o yz = 0. To calculate 
a xy , we cannot directly apply Eq.(31) and Eq.(33) on 
H(k x ,k y ;k z ) for a generic k z , because C4 or Cq rota- 
tional symmetry is not a symmetry of the underlying 2D 



subsystem. However, the subsystems with k z = or 
k z = 7r are invariant under M xy and, on these planes, 
C4 or Cq is restored. The restored rotational invariancc 
allows one to calculate the Chern number on k z = or 
k z = 7r. The insulating condition requires that the Chern 
number on every k z must be the same and the total Hall 
conductivity, a xy , is still given by Eq.(77). In fact, the 
insulating condition requires 



n 



= 1, 



Ci(T)Ci(X)Ci(Y)Ci(S) 
Q{Z)aU)UT)UR) 

&(r)&(M)CipQ _ 



iGocc. 



&(Z)&(A)&(fl) 



11 r}i(A)9i(H)Ci(L) ' 



(78) 
(79) 
(80) 



for n = 2,4,6 respectively 34 . Other values of the above 
products signal the presence of Weyl fermions. 

To summarize our findings concerning the Hall con- 
ductance in 3D insulators with point group symmetries, 
we note that only C n for n — 2,3,4, and 6 , C n h (n = 
2, 3, 4, 6) and S n (n = 2, 4, 6) can have nonzero Hall con- 
ductance quantized to a multiple of the number of layers 
along the third direction. 3D insulators with any other 
point group symmetry must have a xy = a yz = a xz = 0. 
In particular, with the exception of 6*2, or 3D inversion, 
which is compatible with nonzero a xz and c^, all point 
group symmetric topological insulators have vanishing 
a xz and o yz . assuming z-axis to be the principle axis. 



B. Magnetoelectric Effect 

It is known that topological insulators with time- 
reversal or space inversion symmetry in 3D posses a coef- 
ficient of the magnetoelectric term E • B, P3, which must 
be quantized to either zero or one half 32,357 . We now 
desire to find if P3 is quantized in the presence of a gen- 
eral PGS in a 3D insulator. We begin our discussion by 
considering the magnetoelectric effect in a 3D topolog- 
ical insulator with an arbitrary rotation symmetry, the 
operator of which is denoted by R. In 3D, a rotation 
R can either be a proper or an improper rotation. In a 
proper rotation, the system is rotated about a given axis 
in space by a given angle; while in an improper rotation is 
the combination of a proper rotation and space inversion, 
i.e., (x, y, z) — > (—x,—y,—z). Mathematically, if Rij is 
the three-by-three rotation matrix, then det(i?) = 1 for 
proper rotations and det(i?) = —1 for improper rota- 
tions. A general operation in a point group is either a 
proper or an improper rotation. 

By Eq.(12), one can use £>(k) to express \u(k )} in 
terms of R\u(k)), where k' = Rk is k transformed by R: 



|u (k')> = *&(k)£Mk)>. 



(81) 



b^occ 
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Using Eq. (81), the non-Abelian Berry connection Ai(k) This non-Abelian gauge transformation exactly takes the 
then has the following property: form of the non-Ableian transform of a gauge potential, 

with the exception of the presence of the prefactor Rij. 
(A(k')) ab = -i(u a (k')|%fy|u 6 (k')) (82^he transformation of the field strength, whose proof is 

E, „ , . N , - i + , , ~, , .included in Appendix E, is 

c^eocc F ij (k) = R ii ,R jj ,B(k)F i , j ,(k)B^k). (83) 

= -iRij ^2 ^ac(k)(w c (k)|9 i BJ d (k)|ii d (k)) We are now in a favorable position to understand the 

c,<Zeocc constraints on P3 imposed by the presence of general ro- 

= R l - j (B(k)A J (k)B- 1 (k)) ab - i%(B(k)9 J S(k)) afe tation symmetry R. In fact, one has 

I 



+ P z = J^J d 3 ke l]k Tr{(F l3 (k) - ^Mk)A 3 (k))A k (k)} (84) 
x -L J d 3 ^Tr[(^-(k')-^A(k')A(k'))A(k')] 

— l — I d i ke l]k R ll ,R rr R kk ,Tr[(BT ¥r B^ - 2 -^-(BA^ - iBd v B^){BA^ - iBd r B^))(BA k >B^ - iBd v B% 



where n is a gauge dependent degree of freedom. To proceed, notice that 

£ij k Rii'Rj 3 'R kk ' = det(i?)e l / J v A ,/ = ±£j/j/ fe / . (85) 

If R is a proper rotation plus sign is taken and if R is an 
improper minus sign is taken. Then we have 



n + P 3 = y|- / d 3 ke l]k Tr[(B^B^ - j{BA,B^ - iBd^){BAjB^ - iBd^)){BA k B^ - iBd k B^)] (86) 
= it? / ( l'' h ''^ lrU: -- - \iAA)M + / rf 3 ^r[(S5 i Bt)(B9 J Bt)( S9fcS t)] 



= ±(n + P 3 ) ± ^2 f (PkTrKBdiB^iBdjB^iBdkB*)]. 
If P is improper, we have 



Ps = / d3fce »^ Tr [( i3a " Bt )(' B ^' Bt )(' S ^ Bt )] mod !. ( 87 ) 



which is indeed one half of the winding number of B. 
However, if R is proper, the two Pa's cancel each other 
and the symmetry has no constraint on P3. Additionally, 
it gives 

-L J d?kTr[{Bd i B^){Bd j Et){Bd k B^)} = Q. (88) 

Eq. (88) indicates that although one may define the 
winding number of the sewing matrix as a quantum num- 
ber, that number is always zero. 

Up to this point we have formally derived the state- 
ment that only systems with improper rotation symmetry 



can have quantized P3. In fact there is a simple way to 
understand this simple result. In an insulating medium 
with applied electromagnetic field, any point group sym- 
metry of the medium must be preserved if one rotates the 
applied field together, which means that the axion term 
in the Hamiltonian density P3E • B remains invariant un- 
der some improper symmetry. On the other hand, since 
E-B is a pseudoscaler, under any improper rotation there 
is E B > — E • B. Therefore we have P 3 = — P 3 . From 
this equation one may be tempted to obtain P 3 = 0, but 
is reminded that P3 is only well-defined up to some in- 
teger in a compact space. Therefore P3 = — P3 should 
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be interpreted as P3 = — P3 + integer or 2p; = integer. 
This is the intuitive argument leading to the same quanti- 
zation of P3 in insulators having some improper rotation 
symmetry. 

It is easy to check that every point group that contains 
an improper rotation operation must have any of C s , S2, 
or S4 as its subgroup(s), where C s is a group generated by 
just one mirror plane. Therefore the smallest symmetry 
point groups that have quantized P3 are C s , S2 and S4, 
where it should be noted that £2 is not a subgroup of S4. 
But this does not tell us whether it is possible to have P3 
quantized to the nontrivial value of 1/2 in the presence of 
C s , S2 and S4. In Ref.34 and 35, it was made clear that 
there are insulators having inversion invariance without 
time- reversal invariance that still have nontrivial P3. In 
the following, we will give examples of systems having 
only C s or only £4 that have nontrivial P3. We begin 
with a 3D Hamiltonian 

H(k) = sin(^)r 1 + sin(fc. y )r 2 + sm(k z )T z (89) 

+ M(k)r 0! 

where M (k) = 3 — m — cos(fc x ) — cos(fc y ) — cos(fc y ), r = 
1 ® t z , Ti = a z ® t x , T 2 = 1 ® t v and T 3 = a y ® t x . This 
Hamiltonian has time-reversal symmetry and Oh point 
group symmetry. One can add terms to break it down to 
smaller point groups. To do so, we add a magnetic field 
along z-axis and an electric field along x-axis: 

5U = BT 35 + ET W . (90) 

The magnetic field breaks time reversal symmetry and 
all rotation axes except those about the z-axis, while the 
electric field further breaks rotation symmetry about z- 
axis. The only symmetry remaining is the mirror re- 
flection about xy-plane. On the other hand, P3 must 
remain unchanged in so much as the external fields are 
not so strong as to close the bulk gap, because a mirror 
reflection is still an improper rotation and quantizes P3. 
If we choose the parameter < m < 3, we obtain a 3D 
model that has P3 = 1/2 but no symmetry other than 
a reflection about xy-plane, or point group C s . We can 
also add the term: 

6H 2 (k) = BT 35 + tsm{k x )T 15 -tsm{k y )T 25 (91) 

to our 3D Hamiltonian. In Eq. (91), the first term re- 
moves time-reversal symmetry and all rotation axes ex- 
cept z-axis while the second term breaks both C4 and 
M z separately but preserves their combination S4 = 
C4 * M xy . As before, as long as the added terms are not 
large enough to close the bulk gap, P3 remains quantized 
at n + 1/2 if < m < 3. 

While time-reversal symmetry is not the major topic 
of the paper but it is still interesting to discuss how this 



symmetry can change the previous results. Intuitively, 
since time-reversal operation T, like inversion, sends k 
to k, it is similar to an improper rotation symmetry as 
far as P3 is concerned. Generally, a symmetry operation 
T * R is equivalent to an improper/proper rotation if R is 
a proper/improper rotation. Therefore, P3 is quantized 
when TR is a symmetry of the system for some proper 
R. From the statement is also derived a not very obvi- 
ous result: if the system is not invariant under separate 
inversion (P) or time-reversal (T), but is invariant un- 
der their combined operation (P * T), it does not have 
quantized P3, as P * T is proper. 

V. CONCLUSION 

We study several bulk topological invariants in 2D and 
3D insulators with crystallographic point group symme- 
tries, focusing on finding (i) the constraints placed by 
these symmetries on known topological invariants such 
as the Chern number and (ii) if a PGS gives rise to new 
topological invariants. In 2D, we show that the Chern 
number of a C„ invariant insulator are determined up 
to a multiple of n by by eigenvalues of C m at high- 
symmetry points, where m divides n. In P„-invariant in- 
sulators, we show that the Chern number is constrained 
to be zero, while the electric polarization, or the cen- 
ter of charge position is a new topological invariant, the 
value of which can be determined by eigenvalues of C m 
at high-symmetry points, where m divides n. In 3D, we 
show that only C„, C n h and S n invariant insulators can 
have nonzero anisotropic 3D quantum Hall conductance, 
while insulators with all other point group symmetries 
must have zero Hall conductance in every component. 
We also prove that the magnetoelectric susceptibility of 
point group symmetric topological insulators in 3D is 
quantized to or 1/2, i.e., a Z 2 number, if and only if 
the point group contains at least one improper rotation. 
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Appendix A: Details in Eq.(29) 

In Eq.(29), we showed how the determinant of a Wilson loop enclosing a quarter of the BZ (A = TXMYT) can be 
expressed in terms of the sewing matrices, but omitted several steps in which identity operators are inserted. The full 
details arc given here: 

det(<u<(r)|E/XK-(r))) (Al) 
= det ( (m (r) I Uvx Uxm Umy U y r K (T) ) ) 

= det(( Ui (r) I <7 4 (C 4 - 1 Urx C 4 )C^ 1 (<7 4 fewQ" 1 )C 4 Umy U Yr \uj (T) ) ) 
= det((^(r)|C , 4f/ryC 2 - 1 ;7yMC'4;7A/y;7yr|^(r))) 

= det( J2 (^(r)|c 4 |u (r))K(r)|^ rr | U6 (y))( U6 (y)|c 2 |u c (y))( Uc (y)|^ M |u d (M)> 

a,6,c,eZ,e,/eocc. 

(ud(M) \C 4 \u e (M) ) (u e (M) | tW |«/ (F)) <«/ (y) | E/yr I «j (T)» 

= dct(B C4 (r^ryS^ 1 {Y)U YM B Ci {M)U MY U Yl ) 

= dct(B Cl (T)B c 1 2 (Y)BcAM)) dct{U TY U Y M U MY U YT ) 

= det(B Cl (T)Bc 1 2 (Y)BcAM))- 

In the steps, we notice that we insert operators like ^2 aeocc |u (k))(u a (k) as identity operators, which is allowed if 
and only if C4 is a symmetry of the system and the system is fully gapped. The C4-symmetry guarantees the existence 
of another band of equal energy at 64k while the insulating gap guarantees no mixing from the unoccupied bands. 

Appendix B: Proof by monodromy with C3 and Ce symmetries 

In this Appendix, we prove the Eq.(32, 33) using monodromy argument. 

For C3 invariant systems, choose a loop shown in Fig. 1(c). From C3 symmetry, we have (through Eq.(83)) F xy {k) = 
J~xy(R3^) and therefore the determinant of the loop integral shown in Fig. 1(c) is 



det(W A ) = e l2C%/3 , (Bl) 

T— 1 , 

11 using 03 symmetry we noLice uiai, using rjq.^zoj, l 

CzU„» rsCn 1 = U„i which leads to 



where (U\)ij = (ui(T)\U\\uj(T)). Again using C3 symmetry we notice that, using Eq.(23), Ca l U rK »C a = U TK >, and 



J 3 U K"K^3 — U K'K' 

Ux = UrK"U K » K U KK ,U K , T (B2) 
= u rK »u K <> K C 3 U KK »C 3 C 3 U K " r C 3 . 

After inserting identity operators as done in Eq.(Al), the determinant detU\ simplifies as 

det(Z4) = det(U rK „U K » k B{K)U kk ..B~\k' ')U K „ r B{T)) (B3) 
= dct(B(K)B- 2 (K')B(T)). 

In the derivation we have noticed that K and K are the same point by translational symmetry. Then using 
B 3 = (—l) F Ijy occ xn OC c on Eq.(B3) and combining Eq.(Bl), we have 

e i2Cn/3 = (-i)F 6 .(r)0 i (K)0 i (K'). (B4) 

For Cg we consider the loop shown in Fig. 1(d). The determinant of its loop integral is 

det(Z4) = e iC% / 3 . (B5) 

Notice that C^UtmCq = U rM > and C^UmkC^ 1 = U M > K , and we have 

U\ = UtmUmkUkm'Um't (B6) 
= UtmUmkC^UkmC^ 1 C 6 1 UmfCq. 
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Inserting identity operators, the determinant becomes 

det(W A ) = dct(U rM U MK B C3 {K)U KM Bcl (M)U Mr B Ce (T)) (B7) 
= det(S C6 (r)£ C3 (A)^ 2 1 ), 

where we have used that C3 * Cq = C 2 . Noticing that C 2 = (— 1) F , we obtain 

e iC v/3 = -Q (- 1 f Vl( r)o i (K)Ci{K'). (B8) 

i£occ. 

Appendix C: High-symmetry points in a finite system 

In Sec. 11(B), we have derived the relation between the many-body C„ eigenvalue and the C m (m dividing n) 
eigenvalues at high-symmetry points in the BZ, upon the condition that all these points can be taken in our system. 
In a real system, whether a k-point appears on the grid depends on the total number of unit cells along each direction. 
For example, in a ID system, k = it only appears when there are even number of unit cells. In this Appendix, we 
exhaust the possibilities of the presence/absence of each high-symmetry point in a 2D system. 

On a 2D lattice with periodic boundary, a point on the grid can be generally written as 

k=g^i^b 1 + g^^b 2! (CI) 

where n\ = 1, Ni and n 2 = 1, A 2 . For each high-symmetry point, the corresponding (ni, n 2 ) are given by: 

r : (m,n 2 ) = (1,1), X : (m,n 2 ) = (Ai/2,1), Y : (m,n 2 ) = (1, N 2 /2), M : (m,n a ) = (iVi/2, N 2 /2) for n = 2,4; 

K : (m,n 2 ) = (JVi/3,2JV a /3), K' : (m,n 2 ) = (2Ai/3, A 2 /3) for n = 3,6. 

Therefore, the parity of N\^ 2 determines whether a certain high-symmetry point can exist in the system. For n = 3, 4, 6, 
we also implicitly assume that N\ = N 2 = No, because otherwise the many-body system cannot have the corresponding 
symmetries. 

Disappearance of high-symmetry points changes the results in Sec. 11(B) in two aspects. First one must remove the 
C m eigenvalues at those points from the formulas; second there are more k-points grouped into groups of two, three, 
four and six, respectively, which requires correction of the prefactors in, for example, Eq.(38,45, 47). Below we rewrite 
the formulas for different combinations of parities of Ni and N 2 . 

For n = 2, if N± = even and N 2 = odd 



c 2 \<s> ) = (-i) (F ' 1)w °" (W ' 2) J] Ci(r)CiWI*o>, (C2) 

i£occ 



if N\ = odd and N 2 = even 



if JVi = odd and N 2 = odd 



For n = 3, if N mod 3/0 



For n = 4, if jV = odd 



C 2 |$ ) = (-I) '"' 1 '"""'"- 2 ' Y[ b(r)Ci(y)\$o), (C3) 

i£occ 



C 2 \$ ) = (-l f- 1)N T' N ^ J] Ci(T)|*o>. (C4) 



C«3|<I>o) = (-l)^oao{N-l)/3 -Q fl .(T)|$ ). (C5) 



i£occ 



c 4 |$„) = ( -i)^-i)^(^-i)/4 Y[ 6(r)|$ ). (ce) 



For n = 6, if Ao = even and Ao mod 3^0 

C e \$o) = (-lJ^-D^W*-*)/ 8 ]J 7?? (r)C,:(M)|$o), (C7) 
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if Nq = odd and Nq mod 3 = 

C 6 \%) = (-1)^-1)^(^-3)76 -Q r}i(m{ K)\<i> ), (08) 

iGocc 

if N = odd and Nq mod 3^0, 

C 6 |d>„) = (-i^-UJWAr-D/e JJ r/4(r) . (eg) 

Appendix D: Vanishing components of Hall conductance in presence of certain point group symmetries 

With point group symmetries, some components of Hall conductance are necessarily zero. There are three compo- 
nents of Hall conductance o~ xz , a yz and o~ xy . By definition, we have j a = o~ a bEb- Under rotation R : r a = R a br' b , we 
have 

Rabj'b = a a bRbcE' cl (Dl) 

or 

j' a = {R- l aR) ab E' h . (D2) 

R^aabR = R'^a'b'Rb'b, (D3) 



In terms of quantum operators, this means 



crab = lim 1 6 - / j a (t)j 6 (0)e iwt dt. (D4) 



where 

w=o 2w 

Specially, if R is an n-fold rotation about z-axis, we have 

Rr^GxzR = cos(27r /n)a xz — sin(27r /n)a yz , (D5) 
R~ a yz R = cos(2ir/n)a yz + sm(2ir/n)a xz . 

Now we consider an insulating system with n-fold symmetry, and we further assume that the ground state is non- 
degenerate, we have -R|$o) = A(_R) | } , because a non-degenerate ground state must be a ID representation of the 
symmetry It can be easily proved that o~ xz = cr yz = 0: 

o xz = (*o|o-x«|*o) (D6) 

n *— 4 

i=0,...,n— 1 

= - > (cos( )a xz -sm{ )a yz ) 

i=0,...,n— 1 

= 0, 

in which we have used X* r \b, = 1 and the mathematical identity J2i=o n-i C0S ^7T = J2i=o n-i sm ^T = 0- 
Specially, we have that in a 3D system with more than one rotation axis, o~ xz = o~ yz = o~ xy = 0; and in a 2D system 
with mirror symmetry, a xy = 0, because the mirror symmetry can be seen as a two-fold rotation about an in-plane 
axis. 

Mark that this proof only uses the definition of Hall conductance and that the many-body ground state is a singlet. 
The conclusion applies to any interacting system with a singlet ground state. 

For n = even, there exists a simpler proof utilizing the fact that any D ra -symmetric 2D insulator must have zero 
Chern number. If the system is invariant under an even-fold rotation about z-axis, it must also be invariant under a 
twofold rotation about z-axis. Therefore the 2D plane in k-space defined by k x — is a 2D insulator with at least 
Di-symmetry, therefore a yz (k x = 0) = 0. Then since the quantized o~ yz (k x ) must be smooth in an insulator, it must 
be a constant. Hence cr yz (k x ) = for each k x and finally, the total a yz = 0. 
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Appendix E: Transform of field strength T under point group symmetry R (proving Eq.(83)) 



^•(k) = diAjik) - ^A(k) + i[A(k),^-(k)], (El) 

and we look at the Abelian and non-Abelian terms separately, using Eq.(82): 

diAjQt!) - djAiQt!) (E2) 
= R U 'R 3 y[de{BA r B^ - iBd r B ] ) - dy{BA l -B ] - iBd v B^)} 

= R lt ,R jr B{d t ,A r - dfAi^ + /.'.,•//„. \{d v B)A r B^ + BAj>{d v B^) - (d J ,B)AB^ - BAi>d r B^\ 
-iRii'Rjj'idi'BdfBi - d y Bd v B^), 

[A(k'), Aj(k')} = Ra-Rjj, [BAi'B^ - iBd v B\ BAyEfi - iBd y B^\ (E3) 
= H,,-l!,rB A,-.Ar B~ - iR w R jj ,(BAi>d J ,B i + {d r B)A v B ] - Hi,,- h',/U),-BArB + BA 3 >d t >B r ) 
+R ii ,R jr {d i .Bd jl B ] - dfBdi'B^). 

From these two equations, Eq.(83) directly follows. (In the above expressions the argument of k is suppressed while 
that of k' is kept explicit.) 



Appendix F: Basic properties of the sewing matrix 

Every point group symmetry operation has a certain order, i.e., there exists integer n(i?) for which R n ( R } = E or 
R n(R) = E t f where E and E , 

are identity operation and a 2-k rotation, respectively. The distinction between E and E' 
is necessary as we are interested in the representation in Hilbcrt space of a single fcrmion. For fermions with integer 
spins, the representation of both E and E' are J, but for fermions with half-integer spins, the representation for E 
and E 1 are / and — / respectively. For a sewing matrix associated with a point group symmetry R with order n, we 
have 

( Y[ B(J2»k)) tf = <«i(£" k )|£K-i(-R n " lk )> ( u ^-i (■R n " 1 k)|fi|«i„_ a (-R n " a k))...<«i a (i2k)|A|u i (©3l) 

s— 0, .. . ,n— 1 i± ,i2 ,i n Gocc 

= ( Ul (R n k)\R n \ Ul (k)). 
If R n = E, then we have 

II B(R"k) =1n occ .xn occ .- (F2) 

s=0,...,n-l 

If R n = E, we have 

JJ B(R s k) = (-l) F X NocMcc , (F3) 

s=0,...,n-l 

Now we show that B(k) is also a unitary matrix: 

(B(k)Bt( k ))y = £ (u^Rkmu^k^iuA^MRk)) (F4) 

i' (Eocc 

= <u<(iZk)|u,-(iZk)) 
= 

In this equation we have used that the point group operator R is unitary. This is because all symmetry operators are 
cither unitary or antiunitary (Wigner's theorem) and a point group operation does not include either time-reversal 
or charge conjugation. 

Using the sewing matrix, we can represent \ui(Rk)) in terms of |itj(k)) and the sewing matrix. To see this: first 
notice 

R- 1 \u l (Rk))= ]T |« J -(k))<u i (k)|^- 1 |u i (iZk))) > (F5) 

j £ occ 



19 

then multiply R on both sides and obtain 

Mi?k)> = Btj (k)l«i(k)>- (F6) 

Finally, we consider a sewing matrix that is associated with an antiunitary operator R . First we show that the 
sewing matrix is still unitary: 

(B(k)Bt(k)) y = >T ( Ml (i?'k)| J R> rn (k))((^(i?'k)|i?> tm (k)))*. (F7) 

Then use the antiunitarity of we have 

((^(fl'kJI-R'u^k)))* = (u m (k)\R rl \u 3 (R'k)). (F8) 

And from this we obtain S(k)B'(k) = In occ xn occ - 

And for an antiunitary operator, we can also express |uj(ii'k)) in terms of |iti(k)) and the sewing matrix: first 
notice 

R rl \ Ul (R'k)) = ]T |u,-(k))<u,(k)|#~V(#k))) 5 (F9) 

jGocc 

then act 7?' on both sides and have 



i(Rk)) = ]T) ((« 3 -(k)|fl' V'^-R'k)))*^'!^^)) (F10) 
= ]T (^(iZ'lc) l^'lttj (k)> |uj (k)> 

j'Gqcc 

= ]T ^>,(k)). 
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